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ABSTRACT 

We perforin a reconstruction of the cosmological large scale flows in the nearby Uni- 
verse using two complementary observational sets. The first, the SFI++ sample of 
Tully-Fisher (TF) measurements of galaxies, provides a direct probe of the flows. The 
second, the whole sky distribution of galaxies in the 2MASS redshift survey (2MRS), 
yields a prediction of the flows given the cosmological density parameter, fi, and a 
biasing relation between mass and galaxies. We aim at an unbiased comparison be- 
tween the peculiar velocity fields extracted from the two data sets and its implication 
on the cosmological parameters and the biasing relation. We expand the fields in a 
set of orthonormal basis functions, each representing a plausible realization of a cos- 
mological velocity field smoothed in such a way as to give a nearly constant error on 
the derived SFI++ velocities. The statistical analysis is done on the coefficients of 
the modal expansion of the fields by means of the basis functions. Our analysis com- 
pletely avoids the strong error covariance in the smoothed TF velocities by the use 
of orthonormal basis functions and employs elaborate mock data sets to extensively 
calibrate the errors in 2MRS predicted velocities. We relate the 2MRS galaxy distri- 
bution to the mass density field by a linear bias factor, b, and include a luminosity 
dependent, oc L a , galaxy weighting. We assess the agreement between the fields as a 
function of a and j3 = f(Q)/b, where / is the growth factor of linear perturbations. 
The agreement is excellent with a reasonable x 2 per degree of freedom. For a = , 
we derive 0.28 < j3 < 0.37 and 0.24 < /3 < 0.43, respectively, at the 68.3% and 95.4% 
confidence levels (CLs). For /3 = 0.33, we get a < 0.25 and a < 0.5, respectively, at the 
68.3% and 95.4% CLs. We set a constraint on the fluctuation normalization, finding 
erg = 0.66 ± 0.10, which is only 1.5c deviant from WMAP results. It is remarkable 
that as determined from this local cosmological test is close to the value derived from 
the CMB, an indication of the precision of the standard model. 

Key words: Cosmology: large-scale structure of the Universe, dark matter, cosmo- 
logical parameters 
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1 INTRODUCTION 

For 15 years, the problem of large-scale flows of galaxies 
has seen little attention relative to other probes of the large 
scale structure in the Universe. The data on peculiar veloci- 
ties has been difficu lt to obtain, and the results had contra - 
dictory conclusions (|Strauss fc Willicklll995l ; IZaroubill2002l ). 
They are limited to small redshifts (~ 100/i -1 Mpc) at which 
distance indicators can reliably be used. These earlier forays 
into the subject led to disagreements that few people wanted 
to sift through. But in the interval, the data has improved 
dramatically, thus stirring recent activity in the subject. 

Peculiar velocities are unique in that they pro- 
vide explicit information on the three dimensional mass 
distribution, and measure mass on scales of 20 — 
50/i -1 Mpc, a scale untouched by alternative meth- 
ods. Local peculiar velocity data is, in principle, af- 
fluent in cosmological information. Power spectra and 
correlation functions could be derived from the data 
by direct cal culation or by ma xi mum likelihood tech- 
mqucs (e.g. iGorski et al.l [l989l; iJaffe fc Kaiserl fl995l; 
Freudling et all 1 19991; lJuszkiewicz et all I2OO0I ; iBridle et all 
200ll ; lAbate fe Erdogduf feoog). Direct calculation of low or- 
der moments of the flow, such as the bulk motion and the 
shear could also analyzed within the fr amework of cosmolog- 
ical models (e.g. iFeldman et al.ll2Qloh . Some authors claim 
the bulk flow of mixed catalogs of galaxies a rgues there are 
probl ems with ACDM (IWatkins et al.| [2009). but recent re- 
sults (Nu sser fc Davis! 120111 ) show that the SFI++ catalog 
by itself has a large scale bulk flow that is consistent with 
ACDM, and this analysis has smaller error bars. Other, 
perhaps more ambitious, applications could involve an as- 
sessment of the statistical nature of the initial cosmolog- 
ical lar ge scale fluctuations, i .e. whether gaussian or oth- 
erwise (INusser fc Dekellll993l ; iBernardeau et al.lll995T ). All 
these analyses could be performed with peculiar velocity 
measurements alone. 

Here we will be concerned with a comparison of the ob- 
served peculiar velocities on the one hand and the velocities 
derived from the fluctuations in the galaxy distribution on 
the other. The basic physical principle behind this compar- 
ison is simple. The large scale flows are almost certainly the 
result of the process of gravitational instability with over- 
dense regions attracting material, and underdense regions 
repelling material. Initial conditions in the early universe 
might have been somewhat chaotic, so that the original pe- 
culiar velocity field (i.e. deviations from Hubble flow) was 
uncorrelated with the mass distribution, or even contained 
vorticity. But those components of the velocity field which 
are not coherent with the density fluctuations will adiabati- 
cally decay as the Universe expands, and so at late times one 
expects the velocity field to be aligned with the grav ity field, 
at least in the limit of small amplitude fluctuations l|Peeblesl 



ll980l ; lNusser" et al .1 [T99 lh . In the linear regime, this relation 
implies a simple proportionality between the gravity field g 
and the velocity field v g , namely v g oc g t where the only 
possible time t is the Hubble time. The exact expression de- 
pends on the mean cosm ological density parameter SI and is 
given by (|Peebleslll980l ), 



, x 2/(17) , , 



(1) 



Given complete knowledge of the mass fluctuation field S p (r) 
over all space, the gravity field g(r) is 



v)=Gp / dV<5 p (r')^- 



(2) 



where p is the mean mass density of the Universe. If the 
galaxy distribution at least approximately traces the mass 
on large scale, with linear bias b between the galaxy fluctua- 
tions 5g and the mass fluctuations (i.e. S g = b8 p ), then from 
JTJ and @ we have 



H p 
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E 
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where n is the true mean galaxy density in the sample, 
P = f (Q,)/b with / w Q - 55 the linear growth factor l|Linderl 
2005), and where we have replaced the integral over space 
with a sum over the galaxies in a catalog, with radial selec- 
tion function cf)(r The second term is for the uniform com- 
ponent of the galaxy distribution and would exactly cancel 
the first term in the absence of clustering within the survey 
volume. Note that the result is insensitive to the value of 
Ho, as the right hand side has units of velocity. We shall 
henceforth quote all distances in units of kms -1 . The sum 
in equation <[3j is to be computed in real space, whereas 
the galaxy catalog exists in redshift space. As we shall see 
in H3.ll the modified equation, which includes redshift dis- 
tortions, maintains a dependence on Q and b through the 
parameter /3. Therefore, a comparison of the measured ve- 
locities of galaxies to the predicted velocities, v g (r), gives 
us measure of /3. Further, a detailed comparison of the flow 
patterns addresses fundamental questions regarding the way 
galaxies trace mass on large scales and the validity of grav- 
itational instability theory. 

In this paper we shall make this comparison using the 
best presently available data for both the velocity and grav- 
ity fields. The direct comparison of the peculiar velocities 
is fraught with difficulty. Distances to individual galaxies 
are typically uncertain at the 20% level and are furthermore 
subject to considerable Malmquist bias. We shall elaborate 



1 <j>( r ) i s defined as the fraction of the luminosity distribution 
function observa ble at distance r for a given flux limit; see (e.g. 
lYahil et af1ll99ll) . 
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a me thod that was first presented 15 years ago (|Davis et al.l 
1 19961 ) (hereafter DNW96) and which was designed to alle- 
viate most of the observational biases. But the peculiar ve- 



those of others (e.g. Hudson et al.lll995l; Kaiser et al. 1991; 


Yahill Il988l: 


Strauss & Davis! 19881: Branchini et al.ll2001a: 


Nusser et al. 


200ll:IZaroubi et al.ll2002T) were all meant to be 



preliminary and none of their conclusions were compelling. 

Recently ideal data sets have been assembled, thus al- 
lowing a new, definitive analysis of large scale flows. The new 
gravity field is very well described by the nearly whole sky 
Two Mass Redshift Survey (2MRS) (Huchra et al. 2005b), 
and the new peculiar ve l ocity catalog is the SFI+- 1- sam- 
ple i|Springob et al.ll2007l . [20091 ). The 2MRS has previously 
been used to address the g ravity field in conside rable detail 
(e.g. lErdogdu et al.l [200^ ; iFeldman et al.ll2008l ). and some 
effort has gone into the comparison of the 2MRS predicted 
veloci ties versus the SFI++ measured veloc i ties, i n particu- 
lar bv lPike fc Hudson! (120051 ) : lLavaux et all (|2010l ). To date 
nobody has included a proper treatment of the correlated 
noise in the analysis. Here we shall compare the observed 
versus predicted radial velocities, taking into account a full 
error analysis based on a suite of elaborate mock catalogs 
designed to match the 2MRS and SFIH — h data sets. Wc 
shall use a r efinement of the metho d of orthogonal mode 
expansion by INusser fc Davis! l| 19941 ) (hereafter ND94) and 
iNusser fc Davisl l|l995h (hereafter ND95). Analysis of pecu- 
liar velocity data is inevitably plagued by systematics, ran- 
dom measurement errors and sparseness of the data. The 
methods employed here are specifically designed to minimize 
these biasing, thus achieving a robust unbiased comparison 
between the measured SFIH — h and the predicted 2MRS ve- 
locities 

In Sfjwe introduce the 2 MRS and the SFI++ data sets 
and various trims that we do to ensure unambiguous reli- 
able results. In fj3] we describe the method for extracting 
large scale peculiar motions from both data sets. We discuss 
the linear equation for predicting the peculiar velocity field 
associated with a distribution of galaxies in redshift space 
and review our old method deriving estimates of galaxy pe- 
culiar velocities from the inverse Tully-Fisher (ITF) relation 
by means of an expansion over orthonormal modes (basis 
functions). We focus on the new refinements designed to op- 
timize the extraction of the signal from the data. As has 
been the case in the past, mock catalogs constructed from 
iV-body simulations are essential for debugging and calibrat- 
ing the methods. This is especially so for our application, 
since the entire analysis is performed in essentially pure red- 
shift space. We present details of the mocks in Sj4] In Sj5]we 
inspect the flow fields reconstructed from the 2MRS and 
SF+- 1- data, visually and statistically, demonstrating that 
differences between them are similar to those expected in 



the mocks. We present our constraints on the cosmological 
parameters in Sj6] In the concluding iJ7] we summarize our 
findings, discussing their implications and contrasting them 
with other results in the literature. For readers wanting to 
avoid the 'how to' details, we suggest skipping [J3]and [|4]but 
then coming back to understand how our machinery oper- 
ates. 



2 NEW DATA FOR THE COMPARISON 
2.1 Gravity field 

Twenty years ago the only catalog of galaxy photometry 
with uniform c overage over the full sky was derive d from the 
IRAS satellite |Strauss fc DaviJl98Sl ; lYahi]||l988l ). From the 
point source catalog (galaxies were unresolved in IRAS) a 
flux limited sample at 60/im was constructed and redshifts 
were obtained for all objects to construct the IRAS PSCz 
(Point Source Catalog redshift survey, Strauss et al. 1990). 
Among other problems, this PSCz catalog gave little weight 
to ellipticals (which are dim at 60pm as this wavelength is 
dominated by dusty star formation) and suffered from severe 
confusion in regions of high density. However the uniform 
full-sky coverage was unique in enabling the estimation of 
local gravity, and furthermore our local gravity field (in a 
relatively low density region of the universe) is dominated by 
spiral, not elliptical, galaxies and IRAS gave a fair, although 
noisy, representation of the spirals. 

Much larger redsh ift surveys do now exist e.g. SPSS 
l|Abazaiian et al.l 120091 ) and 2dF (jColless et al.ll2003l ), but 
few have attempted to be complete over the whole sky as 
many cosmological measurements do not require such com- 
plete surveys and a trade off has been made between depth 
and sky coverage from the available telescope time and 
resources. The most recent *all sky* ima ging survey was 
the T wo Micron All-Sky Survey (2MASS, ISkrutskie et all 
(2006)), and the 2MASS Extended Source Catalog (XSC, 
Jarrett et al. 2000) extracts from that imaging a flux limited 
(to K=13.5) sample of half a million extragalactic objects. 
The 2MASS Redshift Survey (2MRS, Huchra et al. 2005a) is 
a program to obtain redshifts for all galaxies in the 2MASS 
XSC to a fixed flux limit in the K-band. The K=11.25 mag- 
nitude limited version of 2MRS consists of 23,000 galaxy 
redshifts with uniform sky coverage to within 5 degrees of 
the Galactic plane except towards the Galactic centre where 
stellar confusion limits the catalog to +/-10 deg (Huchra et 
al. 2005b). A K=11.75 mag limited version of 2MRS is al- 
most complete, consists of 43,000 galaxy redshifts, and will 
be made available soon (Huchra et al. in prep.). Since the 
sample is K band selected, the extinction correction is mod- 
est and it is ideal for calculating local gravity. 

In the Southern hemisphere redshifts for the 2MASS 
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galax i es we re observed as part of the 6dFGS (jjones et al.l 
120051 , l2009h . which used the 6dF multi-fiber spectrograph 
on the 1.2m UK Schmidt in Siding Spring, Australia. Their 
ultimate product was a map of 110,256 2MASS galaxies in 
the Southern sky to a magnitude limit of K=12.75 mag 
and to within 10 deg of the Galactic plane. This survey 
is far deeper than the stated goal of 2MRS, but also has a 
higher Galactic latitude limit. In the Northern hemisphere 
the 2MRS builds on a strong tradition of reds hift surveys 
at the CfA: the CfA redshift su rvey and ZCAT (|Davis et al.1 
Il982l : Ide Lapparent et alJll986h . In the absence of a north- 
ern hemisphere equivalent to the 6dF, new redshift obser- 
vations are done galaxy by galaxy using the 1.2m telescope 
at the Fred Lawrence Whipple Observatory on Mt. Hop- 
kins, AZ. The average density of galaxies at the magnitudes 
2MRS is observing is about 1 per degree, so without a wide- 
field multi-object spectrograph in the Northern hemisphere 
this remains the most efficient way to get new redshifts. 
Lower Galactic latitude galaxies in the Southern hemisphere 
(|6| > 5) were added to 2MRS from observations at CTIO. 

The version of the 2MRS which is complete to K=11.25 
(consisting of 23,200 galaxies; Huchra et al. 2005, Westover 
2007) has been used in to calculate t he acceleration on the 
Local Group by (jErdogdu et alJl006al V The dipole estimate 
seems to converge to the CMB result within 60/t -1 Mpc, 
suggesting that the bulk of the motion of the Local Group 
comes from structures within that distance. They also have 
done a dipole analysis, weighting the sample by its lumi- 
nosity, rather than the counts, and find relatively minor 
cha nges. Density and v elocity fields have been calculated 
by l|Erdogdu et alj|2006h for the K=11.25 sample. All major 
local superclusters and voids are successfully identified, and 
backside infall onto the "Great Attractor" region (at 50/i _1 
Mpc) is detected. 

The 2MRS catalog appears to be a fair tracer of the 
underlying mass distribution. The real-space correlation 
lengths, ro is best fit by a regression ro = (7.5 ± 0.5) — 
(3.0 ± 0.6)log in n , where n is the cumulative nu mber density 
in 10 ~ 3 h 3 Mpc (|Westoverll2007f l. In contrast, jZehavi et al.1 
|2010T ) report that the R band optically selected SDSS survey 
gives rp = (6.7±0.1) — (2.0±0.1)log ln rt. In terms of bias esti- 
mates, iWestoyej2007f ) reports 6/6* = 0.73+0.24L/L, while 
(|Norberg et al. II2OO2J) state that 6/6, = 0.85 + 0.15L/L, for 
the 2DF survey. In other words, Westover's data show the 
2MRS correlations are more dependent on luminosity than 
are optically selected samples. In view of this luminosity de- 
pendent result, it makes the most sense to evaluate the grav- 
ity field in a luminosity-weighted manner; i t is computed be - 
low with a variety of luminosity weightings. IWestoverl (|2007l ) 
has also made a mock catalog for the missing galaxies at low 
latitude by interpolating the galaxy density above and be- 



low them in three dimensions. We shall use this catalog as 
an estimate of the local mass density. 



2.2 TF sample 

Twenty years ago, the mis-calibration of full sky Tully- 
Fisher data was the problem that led to very discrepant 
results for the det ermination of ft = Q/6, with fi = 0.5 ± 0. 2 
|Davis et alJll99d ) and p = 1.0±0.2 (e.g. lDekel et~al]|l993l ). 
The mistaken TF calibration led to a large scale flow that 
confused both analyses, but in the end, it was a calibration 
error in the Souther n sky which made a false large-scale flow 
jWillick et al.lll997t ). In one analysis this led to a higher \ 
than was acceptable, and in the other it led to a biased 
result. 

For the analysis below, we use the recently completed 
survey o f spiral galaxies with I -band Tully-Fisher di s tance s, 
SFI++ ^Masters et all |2006| ; ISpringob et~ai1 120071 . l200Sh . 
which in turn builds on the origin a l Spiral Field I-band Sur- 
vey jGiovanelli et al.1 Il994l . Il995l; lHavnes et al.lll999h and 
Spiral Cluster I-band Survev ifGiovanelli et al.ll997al lb). We 
use the published SFI++ magnitudes and velocity widths, 
and derive our own peculiar velocities, rather than taking 
the published distances as given. We use the SFlH — h catalog 
as it includes several datasets to give full sky coverage. It 
is not essential for our analysis that the peculiar velocity 
sample have uniform sky coverage, but they must have a 
uniform calibration. 

The other ma jor Tully-Fisher catalog was published by 
|Tu11v et alj|2008j ). This survey is restricted to cz < 3000 
km/s, and includes many of the same galaxies in that red- 
shift range as SFI++. Tully also make use of a different 
algorithm for measuring spectral line widths, which are not 
easily comparable to the values derived for SFI++. So while 
one could in principle combine the two catalogs for this anal- 
ysis, the small potential gain in sample size is not enough to 
justify the resulting heterogeneity in observational methods 
and data analysis. 

We shall use the inferred distances, as well as red- 
shifts, to derive an estimate of peculiar velocity for each 
galaxy. Correlation analysis (|Borgani et al.l 120001 ) indicates 
that peculiar velocities in the SFI-(— (- behave as ex- 
pected f or ACDM mo d els, ( see esp. (Branchini et al.l 
2001 bl; iFreudling et al.l Il999l ; Ida Costa et al. Il998l ; 
Feldman fc Watkinsll20cia » 



In the analysis below, we shall use the inverse of the 
Tully-Fisher (ITF) relationship, as given in equation [H] We 
begin by drawing t he published magnitudes, ve locity widths, 
and redshifts from lSpringob et al.l (|2007l . 120091 ). We include 
all field, group, and cluster galaxies, which leaves us with an 
initial sample of 4859 galaxies. Galaxies in groups and clus- 
ters are treated as individual objects, though the redshifts 
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Figure 1. top: The scatter of the r\ — Mj relationship, where 
7 = 0.12 Bottom: The mean Ay in .5 magnitude bins of the raw 
distribution (recall that r\ = log(W)). The red is for galaxies with 
cz > 5000 km/s, while the blue is for cz < 5000 km/s. Note the 
change where Mj > —20, and also how the red and blue curves 
appear to have different TF relations. The histograms are for the 
raw data, with no flow model, while the points, with ltr error 
bars, are the mean values of rj — Mj after the fit. The nearby and 
more distant galaxies now have identical TF relations. 



for template cluster galaxies are replaced by the systematic 
redshift of the cluster. Following Giovanelli et al. (1997a), 
we brighten the magnitudes of Sb galaxies by 0.10 magni- 
tudes and brighten the magnitudes of spirals earlier than 
Sb by 0.32 magnitudes, while leaving types later than Sb 
unchanged. This is done in order to account for subtle dif- 
ferences in the TF relation of different spiral subclasses. We 
select only objects with inclination i > 45 deg to ease prob- 
lems with inclination corrections. The data must be trans- 
formed to the LG frame, and galaxies with cz < 200 km/s 
are deleted. All the analysis is done in the LG frame as 
the boundary conditions then simply becomes v g — > and 
v it f -> 0. 

The few galaxies with large residual Arj (A77 is the resid- 
ual from equation [S] see figure [TJ are sufficiently deviant to 
be a worry for statistics which depend on data with a com- 
pact core and no long tails. The typical outlying object is 
not unusually nearby in redshift, and peculiar motions can- 
not be the explanation. The vast majority of SFI++ galaxies 
have well behaved TF relationships; perhaps the outliers are 
undergoing a merger? There are large negative Arj outliers, 
but few corresponding large positive Arj outliers, and this is 
resolved by clipping the outliers at jA^j > 0.20. 



The top panel of figure [T] shows the distribution of Arj 
before galaxies have had their magnitude changed because of 
peculiar velocities; in the bottom panel is shown the result of 
averaging the data into .5 mag bins, where the red histogram 
is for galaxies at cz > 5000 km/s, and the blue histogram is 
for those with cz < 5000 km/s, in the raw data. The red and 
blue points, plus 1 — a error bars, are the average < A77 > 
versus M after the best flow model, described in £13.51 is fit 
to the data. The blue point with M < —23 is deviant, but it 
only represents 11 galaxies, compared to an average of 170 
galaxies in the other bins. Note that the zero point of the 
< Arj > behavior makes no difference; only the slope, the 
constancy of < A77 > versus M is important. Compared to 
the different slopes before the flow corrections are applied, 
the TF relation is now identical in the foreground versus the 
background of the SFlH — h data. This figure is for illustration 
only, as the data is not binned during the fitting process. 

The bend in the TF relationship at M = —20 is known 
to be a result of the reduced mass in the Baryonic T ully- 
Fisher relation (|Stark et al.ll2009l ; iGurovich et alJl201Ch . We 
are missing the data to straighten out the curve, and since 
the 1TF method is easiest to apply if there is a linear rela- 
tionship between rj and M, we simply delete all galaxies with 
AI > —20 from further consideration. After all these cuts, we 
are left with 2830 spiral galaxies with 200 < czlg < 10000 
km/s. 

The raw distribution of Arj, after limiting the sample 
and fitting the best linear curve, is found to approximately 
fit a gaussian with a — 0.059. This is the dispersion with no 
flow model applied. The gaussian width to the distribution 
is a = 0.0558 after the flow model is applied. The small de- 
crease is limited by the intrinsic, dominant noise of the TF 
relation. This noise has numerous causes, such as the uncer- 
tainty in the inclination correction of the SFI++ galaxies, 
or small variations in the outer limits of the rotation curves 
of the galaxies. 

The following sections, Sj3]and iJS] explain the machinery 
for effecting this reduction. 



3 RECONSTRUCTION OF PECULIAR 
VELOCITIES 

In this section we outline our method described in ND94, 
ND95 and DNW96 for deriving the smooth peculiar veloci- 
ties of galaxies from an observed distribution of galaxies in 
redshift space and, independently, from a sample of spiral 
galaxies with measured circular velocities rj and apparent 
magnitudes m. 
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3.1 Peculiar Velocities from the Distribution of 
Galaxies in Redshift Space 



from 
19951) 



There are several methods for generating pecul iar velocities 

redshift surveys, using linear (e.g. Fisher et al.l 

and non-li near rela tions (e.g. iPeeblesl 198C& 

Croft fc Gaztanaeal 1 19981: 



Nusser fc Branchinil 



Frisch et all |2002eJ ; lEnfilin et all 120091 ') Here we restrict 
ourselves to large scales where linear-theory is applica- 
ble. We will use the method of ND94 for reconstructing 
velocities from the 2MRS. This method is is particularly 
convenient, as it is easy to implement, fast, and requires 
no iterations. Most importantly, this redshift space analysis 
closely parallels the ITF estimate described below. We next 
present a very brief summary of the methodology. 

We follow the notation of DNW96. The comoving red- 
shift space coordinate and the comoving peculiar veloc- 
ity relative to the Local Group (LG) are, respectively, 
denoted by s (i.e. s = cz/Ho) and v(s). To first or- 
der, the peculiar veloci t y is i rrotational in redshift space 
l|Chodorowski fc Nusserj Il999l ) and can be expressed as 
V g (s) — — V$(s) where <&(s) is a potential function. As an 
estimate of the fluctuations in the fractional density field 
So(s) traced by the discrete distribution of galaxies in red- 
shift space we consider, 



*o(«) = 



1 \ -> w(Lpi) 



exp 



(s - Si) 



2<r 2 



1 .(4) 



l/4> = 3.7 times the weight of foreground galaxies in the 
generation of the velocity field, v g . 

If we expand the angular dependence of "1? and So(s) 
redshift space in spherical harmonics in the form, 



200cl ; ®(8) = ^2^2 $i m (s)Y lm (9,tp) 



(•>) 



!=0 m= 



-I 



and similarly for So, then, to first order, $; m and <5oim satisfy, 

2 d$i m \ 1 l(l + l)$lm 



-— (< 

s 2 ds \ 



ds 



1 + 



(6) 



(Solm - K(s) 



d$ 



I m 



ds 



where 

dln0 2 dbaw(Loi) 



ds 



,, , (7) 

represents the correction for the bias introduced by the gen- 
eralized Kaiser rocket effect. As emphasized by ND94, the 
solutions to equation ((6J for the monopole (I — 0) and the 
dipole (I = 1) components of the radial peculiar velocity in 
the LG frame are uniquely determined by specifying vanish- 
ing velocity at the origin. That is, the radial velocity field 
at redshift s, when expanded to harmonic I ^ 1, is not in- 
fluenced by material at redshifts greater than s. 

In this paper, we shall consider solutions as a function 
of /3 and the parameter a defining a power law form Wi oc Lf 
for the galaxy weights. 



where n — ^2 . w(Loi) / 4>(si) and w weighs each galaxy ac- 
cording to its estimated luminosity, Loi. The 2MRS density 
field is here smoothed by a gaussian window with a redshift 
independent width, a — 350 kms -1 . This is in contrast to 
DNW96 where the IRAS density was smoothed with a width 
proportional to the mean particle separation. The reason 
for adopting a constant smoothing for 2MRS is its dense 
sampling which is nearly four time higher than IRAS . We 
emphasize that the coordinates s are in observed redshift 
space, expanded in a galactic reference frame. The only cor- 
rections from pure redshift space coordinates is the collapse 
of the fingers of god of the known rich clusters prior to the 
redshift space smoothing (Yahil et al. 1991). Weighting the 
galaxies in equation Q by the selection function and lumi- 
nosities evaluated at their redshifts rather than the actual 
(unknown) distances yields a biased estimate for the density 
field. This bias gives rise to Kaiser's rocket effect (|Kaisej 
1 19871 ). 

To construct the density field, equation [4j we volume 
limit the 2MRS sample to 3000 km /s, so that <j>(s < 3000) = 
1, resulting in tf>(a = 10000) = 0.27 (IWestoverl2007l ). In prac- 
tice, this means we delete galaxies from the 2MRS sample 
fainter than M* +2. Galaxies at 10,000 km/s therefore have 



3.2 Peculiar Velocities from the Inverse 
Tully-Fisher relation 

Given a sample of galaxies with measured circular velocity 
parameters, T)i = logoj^, linewidth un, apparent magnitudes 
rrii , and redshifts Zi , the goal is to derive an estimate for the 
smooth underlying peculiar velocity field. We assume that 
the circular velocity parameter, r), of a galaxy is, up to a 
random scatter, related to its absolute magnitude, M, by 
means of a linear inverse Tully-Fisher (ITF) relation, i.e., 

n = jM + rio. (8) 

One of the main advantages of inverse TF methods is 
that samples selected by magnitude, as most are, will be 
minimally plagued by Mal mquist bias effec t s when analyzed 
in th e inverse direction l|Schechterl Il980l ; lAaronson et al.l 
1982). We write the absolute magnitude of a galaxy, 

Mi = M 0i + Pi (9) 

where 

M 0i =rm + 51og(zi) - 15 (10) 
and 
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Pi = 51og(l - Ui/zi) 



(11) 



where m; is the apparent magnitude of the galaxy, Zi is its 
redshift in units of kms -1 , and m its radial peculiar velocity 
in the LG frame. 

ND95 base a velocity model on spherical harmonics and 
spherical Bessel functions, for galaxies distributed over the 
sky to 6000 km/s. With the 2MRS we extend the gravity 
field to 10,000 km/s. In general, one can write the function 
Pi in terms of an expansion over j m orthogonal basis func- 
tions, Ff, 



E"' ' 



3 = 



with orthonormality conditions, 



E' 



(12) 



(13) 



N g , where N g 



and the zeroth mode defined by F° — 1/ 
is the number of galaxies in the sample. The mode F° de- 
scribes a Hubble-like flow in the space of the data set which 
is degenerate with the zero point of the ITF relation. Here 
we set a — 0, which removes the Hubble-like flow from the 
gravity field, below. The best fit mode coefficients, a\ the 
slope, 7, and the zero point 770, are found by minimizing the 
X 2 statistic 



2 

XiTF 



E 



(7M0; + j Pi + 770 - m) 



(14) 



'?;,int 



where cr^int is the rms of the intrinsic scatter in 77 about the 
ITF relation, and iV g is the number of galaxies in the sam- 
ple. Given the orthonormality condition, the solution to the 
equations dx 2 TF /da j = 0, dx 2 TF /dj = and dx 2 TF /drjo = 
is straightforward. Thanks to the orthonormality condition, 
the covariance matrix < SajSa'j > of the errors in a? is di- 



agonal with 
a a =< (So?) 1 



> 



1/2 . 



7 



(15) 



This lack of covariance of the errors in the coefficients is most 
rewarding as it makes the ITF error analysis exceptionally 
simple. Therefore, statistical assessment of the match be- 
tween the data will be done at the level of the modes rather 
than the peculiar velocities. The interested reader will find 
details in ND95 and DNW96. 

3.3 The orthonormal basis functions 

The choice of radial basis functions for the expansion of the 
modes can be made with considerable latitude. The func- 
tions should obviously be linearly independent, and close to 



orthogonal when integrated over volume. They should be 
smooth and close to a complete set of functions up to a 
given resolution limit. ND95 chose spherical harmonics Y~ ( m 
for the angular wavefunctions and the derivatives of spheri- 
cal Bessel functions for the radial basis functions, motivated 
by the desire to use functions which automatically satisfy 
potential theory boundary conditions at the origin and the 
outer boundary. That is, they chose 



P(y,t 



m— I 

_ E 

n=0 Z = m— — l 



E 



Q>nln 

y 



{i'i{Ky)-c n )Y lm {0,<t,) .(16) 



The function y(z) is designed to compress the distance 
scale, increasing the smoothing scale of the mode to deal 
with increased noise at large distances. For this analysis we 
use 



(z) = v/log(l + (z/Zsf) 



(17) 



where z s — 5000 km/s. The constant c.\\ is non-zero for 
the dipole term only and ensures that P = at the origin, 
and is non-zero at the outer boundary. Details of how the 
orthogonalized functions F? are derived from this expansion 
are given in ND95. 

The spherical harmonics are expanded to a maximum 
n = 5 and 1 = 3, except we delete the n = 1 mode for 
/ = as this mode can be confused with the false Hub- 
ble flow described in the next section. We also include an 
external quadruple, distinct from the internal quadruple, 
to describe the gravity induced by material at distances 
cz > 10000 kms -1 . Summing over the values of m, that 
makes a total of 72 modes fit toward reducing the x 2 of equa- 
tion !14l The use of y(z) is designed to allow the radial resolu- 
tion to degrade with distance; for example, the n — 5 modes 
have a wavelength in the radial direction of 3000 kms -1 
at cz = 8000 kms -1 and a wavelength of 1300 kms -1 at 
cz = 1000 kms" 1 . 



3.4 Expanding the 2MRS gravity field 

In order to assess the match between the velocities by means 
of the expansion coefficients and to ensure that both fields 
are smoothed similarly, the 2MRS predicted velocities v 
must also be described by an expansion over the basis func- 
tions used in the w itf model. Using the machinery for com- 
puting a gravity field described in §3.11 one can generate 
a linear theory predicted peculiar velocity v g for any point 
in space as a function of its redshift for any value of /3. We 
must ensure that the smoothing scales of the ITF and 2MRS 
predicted peculiar velocities are matched to the same reso- 
lution. Therefore we expand v & in terms of the modes used 
in the velocity model. Because of the orthonormality, we can 
write the mode coefficients as 
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(v gi - H'czi) 



(18) 



where the H' term is a correction for Hubble flow and the 
summation index i is restricted to be over the positions of 
the same galaxies in the ITF expansion. 

This procedure will filter out fluctuations that are not 
described by the resolution of our basis functions. We do not 
include any mode such that Pi = constant, which would be 
a pure Hubble flow. In the fitting for the ITF modes, pure 
Hubble flow is absorbed into a shift of the zero point rjo 
and the orthogonality is ensured. Within a given set of test 
points occupying a volume smaller than that used to define 
the gravity field, it is possible for v g to have a non-zero value 
of Hubble flow H' , which must be removed from v g before 
we tabulate the mode coefficients. That is, we tabulate the 
mean Hubble ratio 



H' = 



EN, 
fcfl Vg,kCZ k 



(19) 



and subtract it from the predicted field v g . This "breathing 
mode" which mimics a Hubble flow is not trivial in ampli- 
tude, and can be a 10% correction on the effective Hubble 
constant within simulated catalogs. This mode is cosmolog- 
ically expected to a modest degree, but a bigger portion of 
the effect is caused by error in the determination of n, which 
we estimate by assuming the weighted counts within 12,000 
km/s is the mean value. 

3.5 Refinements of the ND95 functions 



TheE usscr fc Davisl (|l995h reconstruction of the base func- 
tions can only provide a rough estimate of the spatial distri- 
bution of galaxies in the TF data. It does not guarantee that 
the signal to noise in the filtered fields is uniform all over the 
sample. The ND95 expansion yields reliable TF velocities of 
nearby galaxies, but very noisy estimates at larger distances. 
Moreover, it is difficult to achieve a desired resolution as a 
function of redshift and to ensure equal resolutions in the 
radial and angular directions. The ND95 method expands 
the observed velocities in terms of harmonic functions, but 
the individual harmonic modes are not regularized and may 
acquire unrealistically large amplitudes, depending on the 
spatial coverage of the data. 

We aim here at generating base functions which are 
themselves smoothed with a variable isotropic smoothing 
window designed to yield a constant signal to noise in the 
estimated v itl . We construct these new basis functions with 
the help of the ND95 orthogonal functions denoted here by 
F lND . Suppose a single radial velocity field, V aaod , with the 
appropriate variable smoothing has been found. We term 
V , the seed field as the new modes will stem from it. As 



will be described below this field will be chosen as the 2MRS 
predicted velocity field, but any other field representing a 
viable velocity field could serve as V sccd . Given 



5iog(i-y seedi< M) 



(20) 



(where i refers to galaxies in the TF sample) we expand 
P eeed j in base functions constructed according to ND95, 



■.F 



NDj 



(21) 



Here, the number of the modes F ND is sufficiently large so 
that the inverse transformation ^ . o? iaBd F^ D: ' reproduces 
P sood i . In practice we use about 1400 ND95 modes (we go 
to / = 17). 

We then form additional fields, P a , according to 



p . 

± a, i 



NDj 



(22) 



where 1Z is a set of normally distributed random numbers 
with zero mean and standard deviation of unity. This recon- 
struction of the additional fields preserves the "power" in 
the modal expansion and randomizes the phases. So far all 
these fields, P a ,i, are unfiltered and may contain non-linear 
small scale fluctuations. Therefore, we smooth all fields P a ^ 
according to 



= P a ,i'W{s lti ,,R a ^,) , 



(23) 



all galaxies 



where S;y is distance (in redshift space) between the galax- 
ies i and i', and W is a gaussian window of width R St i which 
depends on the galaxy i. The smoothing width R s ^ is tuned 
such that the expected error in the ITF velocity of galaxy i 
is ~ 150 kms -1 and therefore it depends on the redshift and 
the local density of galaxies near i. The smoothing length 
at the positions of galaxies in the SFH — h sample is shown 
in figure @; it varies roughly linearly with redshift, rang- 
ing from lh~ l Mpc nearby up-to 30/i~ Mpc for galaxies at 
redshifts ~ 10000 kms -1 . The new basis functions are then 
obtained by orthonormalizing P™ ooth . We refer to the new 
functions by the standard notation F- , with Fl being the 
smoothed -P cod 4 . These new modes will be used in the ex- 
pansion given in equation (|12[) . In £|3 . 6 1 we describe how we 
determine the number of modes, j m , to be used in the expan- 
sion. The seed field V scod could be constructed by interpolat- 
ing any unfiltered cosmological velocity field on the positions 
of the galaxies in the SFI++ sample. Nevertheless, we could 
improve on this by constructing K. eed from the unfiltered 
2MRS velocities given directly by the solution to equation 
(JBJ. In practice, we use unfiltered v B obtained with /3 = 0.2. 
The choice of u g for /3 = 0.2 is arbitrary; the predicted field 
with any other /3 could be used. If the V r sood ~ v g and v itl 
both describe the same underlying velocity field then the 
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Figure 2. The width of the gaussian smoothing gaussian window 
versus galaxy redshifts in the SFH — h sample. The scatter reflects 
the angular variations in the density of galaxies. 

additional modes should mainly reflect the covariance of the 
errors between the two fields. 

The flow patterns of 9 of the modes are shown, re- 
spectively, in the 9 panels in figure Q. The color scheme 
throughout this paper is normal: red (blue) means outflow- 
ing (infalling), from the central point.. We have extensively 
tested a broad variety of choices for the the first mode. None 
of the results of the analysis reported are sensitive to this 
choice of the first mode. The figure shows that the higher 
order modes exhibit smaller scale structures. This is a direct 
result of the orthogonolization processes. The j th + 1 mode 
has to be orthogonal to the all previous j modes. In order 
to achieve that, the orthogonalized j th + 1 will pick more of 
the small scale structure. 



3.6 The ITF scatter and the number of modes 

Once the basis functions for the modal expansion are given, 
we proceed to solve for the coefficients a? in Pi = J"] a? F- 
by minimizing Xitf m (|14|) . The minimization is also per- 
formed, at the same time, with respect to the slope and 
zero point of the ITF. The estimated slope is 7 = — 0.1297 ± 
0.0015 and -0.13±0.0016, respectively, for 20 and 30 modes 
used in the flow model. The raw slope before fitting the 
model is 7 = —0.1267 ±0.0016. The zero point plays no role 
at all here and we do not keep track of its estimated values. 
All estimated parameters, including a 3 , are independent of 
the (assumed constant) intrinsic scatter a v ,int in the ITF. 
The velocity model can be used to estimate the unknown 
value of o-^jnt. Given the residual 



Ar/i = r)i - (7M 0i + jPi - r)o) (24) 
we approximate o^mt by 

o\ = ^(Arfef/^.o,- (25) 

i 

where 

^.o. t = N B - (j m + 2) (26) 

is the number of degrees of freedom taking into account that 
the minimization of Xitf i s done with respect to j m coeffi- 
cients plus the slope and zero point of the ITF. This a v will 
decreases as the number of modes, j m , in the expansion is in- 
creased. If j m is too large then the higher order modes will be 
dominated by noise. If j m is too small then the model may 
miss significant components of the underlying true galaxy 
velocities. 

The optimal range of j m for our comparison can be seen 
by inspecting the behavior of <j v as a function of j m . The 
(blue) circles in the top panel of figure [5] shows tj v versus the 
number of modes for modes generated from the seed field 
Kecd = V S (P = 0.2) (see 353)) . Most of the reduction in a\ 
is already achieved the first mode. This is very encouraging 
since this means that v s picks up a significant contribution 
of the velocities as described by the ITF data. It also means 
that both the 2MRS and the ITF data are likely to provide 
approximations to the underlying flow field. However, the 
2MRS predicted field v g deviates from the underlying field 
by the presence of correlated errors in the reconstruction 
scheme. The inclusion of additional expansion modes in the 
ITF velocity model will dissolve these errors. The average 
reduction in the variance per mode becomes insignificant 
beyond j m = 64; the average reduction per mode for the 
first 64 modes is 10 times larger than that for the next 75 
modes. An F-test also confirms that the reduction in the 
variance marginal beyond j m = 64. Hence we will approxi- 
mate CT^ int = 0.0558, the value acquired by a v for 64 modes 
in the velocity model. Therefore, the F-test argues that 64 
is the maximum number modes needed to model the ITF. 

The next step is to determine the minimum number 
of modes needed to describe the ITF flow assuming that 
Cj7,int = 0.0558. To do so we tabulate Xitf a function of 
j m and compute the probability Q = Q(xttfI-^, w ) that the 
value Xitf is exceeded by chance (c.f. §6.2 in (Press et al.l 
[l992j). The values of of XiTF/^dot an d Q are represented 
as the (blue) circles, respectively, in the middle and bottom 
panels of figure [5] For j m = 64, we get Xitf /N a o{ = 1 
corresponding to Q — 0.5, in accordance with our choice 
of (7,7, int- Without a velocity model, i.e. j m — 0, we get 
(XiTF/W dof ,Q) = (1-125, 3 x 10" 6 ). This exceedingly low Q 
rejects a vanishing velocity field with very high confidence 
level (CL). Including the first mode alone gives a highly sig- 
nificant improvement: (xitf/^ o c > Q) ~ (1-031,0.12). The 
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Figure 3. The flow patterns of 9 modes for galaxies within 5/i 1 Mpc of the supergalactic plane. The order of the mode is indicated in 
the corresponding panel. 



hypothesis that Xitf value corresponding to the first mode 
is obtained by chance is rejected only at the 0.12 CL. This 
is encouraging since the first mode velocity field is propor- 
tional to the 2MRS predicted velocities, v g (for j3 = 0.2). 
For j m = 30 and 20 we get (x^tf/JV^., , Q) = (1.01,0.29) 
and (1.014,0.2), respectively. 

For comparison with V st , cd = v , the corresponding re- 
sults for random choice of the seed field, V sood , are shown as 
the (red) crosses in the top and bottom panels of figure [5] 
With a random V Beed , about 20 modes are needed to reduce 
a n to the level achieved by the single mode v . 



3.7 What is the purpose of the ITF machinery? 

The expansion of the gravity field is conceptually very clean 
when computed in the LG frame (ND95). The Poisson-like 
equation for the 3-D gravitational field has been solved as 
a sum over the spherical harmonic functions Yi m times 1-D 
functions of r that satisfy physically reasonable boundary 
conditions at the origin. For the purposes of the ITF so- 
lution, we furthermore quantize the radial solutions with 
quantum number n. 

The ITF method is backward from the usual method- 
ology of TF applications; one does not fit curves to the 
scatter of peculiar velocities. Instead, the \ 2 equation [14] 
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Figure 4. Scatter plots of i> itf versus v g (expanded in 20 modes) for galaxies in 9 mock catalogs. About 800 galaxies are plotted for 
each mock. 



is minimized by the addition of linear combinations of the 
orthonormal functions of n,l,m, where each describes a set 
of large-scale flow that satisfy the boundary conditions. Fur- 
thermore, we have endeavored to form a first, 'seed', mode 
based on linear growth rate, but in which large scale grain- 
iness is filtered out by the use of figure [2] 

The individual galaxy's peculiar velocity enters by equa- 
tion [11] with differential dPi oc dui/zi, and since the uncer- 
tainty of peculiar velocity m is proportional to redshift, the 
uncertainty of Pi is redshift independent. This means that 
each object is given equal weight in a fit, and our window 
function is therefore equivalent to the display of figure 
which shows the positions of the SFI++ galaxies. 



4 MOCK CATALOGS AND ERROR ANALYSIS 

As a measure of the agreement between the TF and pre- 
dicted velocities by means of the corresponding expansion 
coefficients, a it( . and a g , we will consider the \ 2 function 



(j,*0=i 

The parameters a and /3 will be obtained by minimization of 
this function. The covariance of the residual, o ltf — a g , is the 
sum of the covariance matrices of the errors in the estimation 
of a itf and a g , respectively. Thanks to the orthonormality of 
the basis functions, the error covariance in the estimation 
of a itt is diagonal with constant terms a\ = (0^/7) 2 (see 
JO] and DNW96). The matrix £ g represents the covariance 
of the errors in the determination of o g . The origin of these 
errors is as follows: 

(i) Equation |S] is expected to predict reliable velocity 
fields only for small amplitude fluctuations. Small scale non- 
linear deviations from linear theory inevitably leak to large 
scales. 

(ii) The 2MRS is a finite number sampling of the under- 
lying density field. This leads to "shot-noise" errors in the 
estimation of the density field. 
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Figure 5. Top: The rms value, cr, ; , of the ITF as a function of 
the number of modes used in the velocity model, middle: Xitf P er 
degrees-of- freedom versus the number of modes. It is unity at 64 
modes. Bottom: The probability that the Xitf excee ds a certain 
value by chance, as a function of the number of modes. The chi- 
squarc is computed assuming an intrinsic scatter iat = 0.0558. 
The value Q = 0.5 is achieved at 64 modes, (see text for details) 
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Figure 6. The peculiar velocities v itl versus t> g of galaxies in 
the real data (top) and the mocks (bottom) for 20 and 70 ex- 
pansion modes, as indicated in the figure. The real v g has been 
reconstructed with = 0.35 and a = 0. Velocities of about 400 
randomly selected galaxies are plotted in each panel where each 
mocks is represented by about 25 galaxies. All velocities in the 
top panel are normalized by the rms value, <r 9 = 233 kms - 1 of 
v , while velocities in the bottom panel are normalized by the 
rms value of v of their corresponding catalogs. 



(iii) Small scale random motions of galaxies, especially in 
groups and clusters, give rise to a smearing of the distribu- 
tion of galaxies along the line of sight in redshift space. 

(iv) There is a possible large scale stochastic biasing 
jDekel fc Lahavl 1 19991 ; ISigad et all |200d : IWild et all 120051 ) 
between the galaxy distribution and the mass fluctuations. 



The only way to achieve a reliable estimate of £ g taking 
into account all of these complicated errors is by means of 
mock catalogs designed to match the general properties of 
the 2MRS. A parent simulated catalog of the whole 2MASS 
catal og has already been prepared ()De Lucia fc Blaizod 
l2007h by incorporating semi-analyt ic galaxy formatio n mod- 
els in the Millennium simulation ( Springel et al.ll2005l ) . From 



this parent catalog we have drawn 15 independent mock 
2MRS catalogues satisfying the following conditions: 

(i) The "observer" in each mock is selected to reside in a 
galaxy with a quiet velocity field within 500 km/s, similar to 
the observed universe. That is, the central server sees only 
one cluster that has high enough peculiar velocities to result 
in negative redshifts. Recall that in the LG frame, the only 
galaxies with negative redshift are in the Virgo cluster. 

(ii) The motion of the central galaxy is 500 to 700 km/s. 

(iii) The density in the environment of the local group, 
averaged over a sphere of 400 km/s radius , is less than 
twice the normal. 
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Corresponding mock ITF catalogs were also prepared. A 
counts-in-cells statistics shows that the distribution of galax- 
ies in the mocks is unbiased relative to the dark matter, i.e. 
6=1. 

The preparation of the mocks for velocity reconstruc- 
tion is done in the same way as the real data. Equation © is 
used to generate prediction of mock v g with /3 = /(f2, A) /b = 
0.47 corresponding to b = 1 and Q = 0.25 and A = 0.75 as 
in the Millennium simulation. The mean of the rms values 
of v g in the mocks is 269 kms -1 and the standard deviation 
from this mean is 56 kms -1 . For contrast, the rms value of 
v g derived from the 2MRS with (3 = 0.35 is 233 kms" 1 . 

To better illustrate the covariance between the residuals 
in the mocks, we plot in figure [4] the velocities v itt versus v g 
for 9 individual mocks. These are velocities expanded with 
20 modes. Note the similarity between the structure of the 
distribution of points in the individual panels and the top 
panel in figure [51 showing w itf versus v (with ft = 0.35) for 
the real data. 



4.1 The error covariance matrix 

The covariance matrix £ B is computed from the 15 mocks 
by projecting the correlation function, £ p , of the residu- 
als AP = Pj t f — P g onto the basis functions, where P lt { = 
51og(l — v it! /cz) and correspondingly for P g . That is 

f g (j, k) =< (<4 - a£)(<& - a k s ) >= £ F^ p (i, i')Fp , (28) 

where the summation in the last term on the right-hand 
side is over all data galaxies and F are the basis functions 
used for the real data. In this calculation, the v itf velocities 
are reconstructed from an ITF relation without adding the 
internal scatter of the TF relation. The reason is that the 
error in a itf resulting from the intrinsic scatter has a simple 
analytic form given by a a - The function £ p is computed from 
the 15 mocks as follows. Denote line-of-sight and projected 
separations in redshift space by, sn and s±, respectively. For 
each mock we tabulate the average < AP1AP2 > over pairs 
with separations defined by the grid. We then normalize this 
quantity by the variance of v it{ (for zero ITF instrinsic scat- 
ter) in the corresponding mock. This is reasonable since the 
rms values of the velocity field vary considerably among the 
mocks and some of them are significantly different than the 
real data. To minimize this cosmic variance and to derive £ g 
given the observed rms value of the velocity, this normaliza- 
tion of < AP1AP2 > for each mock is prudent. The average 
over all mocks is then computed and interpolated from the 
grid onto the actual pair separations in the TF catalog to 
obtain the normalized £ P . The normalized covariance matrix 
£ B is then computed according to (j 281) and scaled by a factor 
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matching the velocity variance estimated from the observed 



5 THE RECONSTRUCTED VELOCITIES 

This section presents a visual inspection of the fields, and 
assess the coherence of the residual v it! — v g by means of a 
velocity correlation analysis. The quantification of the agree- 
ment between the fields and the extraction constraints on a 
and j3 will be deferred to later sections. 

5.1 Visual inspection of the flows 

Blue dots and red crosses of figure [6j show velocities ex- 
panded in 20 and 70 modes, respectively. In the bottom 
panel, v itf versus v g from the mocks are shown. In mock v itf 
velocities are obtained from a fake ITF relation with an in- 
trinsic slope 7 = —0.1 and a scatter with a v — 0.05. For the 
sake of clarity only 400 randomly selected galaxies are shown 
in either panel. Further, each mock is represented by about 
25 galaxies (randomly selected). The velocities are scaled 
by the corresponding rms value of v g in the corresponding 
catalog. 

There is a an excellent overall agreement between v iti 
and v g for 20 modes, both in the real data and the mocks. A 
good agreement prevails even for 70 modes despite the clear 
enhanced noise contamination. In the real data, the rms of 
w itf — v g for 20 expansion modes is 99 kms -1 , significantly 
smaller than a 3 . For 70 modes the rms values of v g and 
w itf — v g are 238 kms -1 and 231 kms -1 , respectively. Both 
panels show clear structures in the distribution of points, 
implying strong covariance between the residuals, v itf — v g , 
in the real, or mock, data. Because the bottom panel repre- 
sents random selections of galaxies from all the mocks, the 
covariance pattern between the velocities is diluted in the 
distribution of points. The covariance pattern is, however, 
clear in figure [3] where scatter velocity plots for a few mocks 
are shown individually. 

In the aitoff projections in Figure [7| we plot the TF 
peculiar velocities, v ttf and the derived gravity modes, v g , 
for galaxies in redshift shells, cz < 2000, 2000 < cz < 4000, 
4000 < cz < 6000, and 6000 < cz < 10000 kms -1 . The 
projections are in galactic coordinates centered on I, b = 
and with b — 90 at the top. Figures [7| and [8] show v with 
P = 0.35. The rightmost plots are the residuals v ;tf — v . 
The key point is to note that the residuals are small for the 
entire sky and have amplitude that is constant with redshift. 
The amplitude and coherence of the residuals v iti — v g is the 
same as for the mock catalogs in figure[Hl where for example 
the lower right picture shows v itf — v g for a mock catalog. 
It is not very dissimilar from the real plot of v itf — v g in 



© 2010 RAS, MNRAS 000, 



14 Davis et al. 




Figure 7. The derived peculiar velocities u tf , v , and u tf — v of galaxies on aitoff projections on the sky in galactic coordinates. The 
rows correspond to galaxies with cz < 2000, 2000 < cz < 4000 4000 < cz < 6000 kms -1 and 6000 < cz < 10000 km/s, respectively. 
The size of the symbols is linearly proportional to the velocity amplitude (see key to the size of the symbols given at the bottom of the 
figure). In order to better see the differences, a 400 km/s dipole, in the direction of the CMB dipole, has been subtracted from the v, tf 
and v velocities. Note that i> itf — v s is considerably smaller than u tf or v g , even for the most distant galaxies. 



the upper right, demonstrating the feasibility of the entire 
method. 

Note the quadrupole pattern for cz < 4000 kms -1 
in figures [7] and [8] visible after 400 kms -1 has been sub- 
tracted from th e flow . This has been previously noted by 
iHaugballe et al l (|2007l ) on the basis of the flows detected in 
133 SNe. The quadrupole is the typical pattern observed in 
nbody simulations and is the principle mode of collapse to 
a 1-D structure. 

There is amazingly overall good agreement between the 
large scale motions as described by v it! and v g . The residual 
velocities are coherent over large scales but they are clearly 
of smaller amplitude than v itt and v g . Note that residuals 



shown in v.. 



particularly visible in the shell 6000 < 



cz < 10, 000 kms , are dominated by I — 4, because the fit 
for the reduction of the TF \ 2 is limited at I — 3 modes. 



5.2 Correlations 

The residuals, both in the real and mock data, have error 
fields, u tf — v , that show large regions of coherence. To ad- 
dress the significance of these errors, we show in figure [9] the 
velocity correlation function, (jGorski et al.|[l98^ ). defined as 



*(s;u) 



E«ltl2COs9l2 
pairs 1 * ^ A 

ECOS 2 #12 
pairs 



(29) 



© 2010 RAS, MNRAS 000, 



Gravity versus Velocity 15 



_^ -6000 
>- 

o 

00 4000 




-6000 -4000 -2000 



2000 4000 6000 



SGX km/s 



Figure 8. Supergalactic plane projection, \SGB\ < 30, of the derived flows. To better see the differences in the plots, a dipole of 400 
km/s towards the CMB pole has been subtracted from the fields, and is shown in the bottom left and bottom center. A quadruple 
velocity is now visible in the plots. The points are drawn at the estimated distance of an SFH — (- galaxy, and the line, blue or red, is 
drawn to the galaxy's redshift. In other words, the length of the arrow is the peculiar velocity. The lower right plot shows u tf — v for a 
mock catalog, and the upper right shows v iti — v g for the data. They have very similar degrees of coherence. 



where the sum is over all pairs, 1 and 2, separated by vec- 
tor distance S12 (in redshift space), #12 is the angle between 
points 1 and 2, and u is either v m (dashed red) or v it! — v 
(red for data, blue for 15 mock catalogs), At small lags for 
the real data, the function *&(r; u tt — v ) is a factor of 3 less 
than ty(s;v it! ), about the same as for the mock catalogs. 
Note how the large coherence of v it! is enormously dimin- 
ished in \[>(s > 2000km/s; w ltf — v ). This shows that the 
coherence seen in the residual field, figure[9] is expected and 
is not a problem. The large scale drift of a sample is demon- 
strated by the persistent amplitude of *I> beyond « 60 — 80 
Mpc. 

The bottom panel of figure [9] shows velocity correlations 
for 15 mock catalogs where the actual velocity, Vtrue, gen- 
erated in the nbody code and then smoothed with the 20 
mode expansion, can be compared to either v itt or v s . Note 
that the raw velocities, vuf (red), have enormous correlation 



that reaches large lag, while the correlations, (vtrue — v ltf ), 
(blue) are extremely small. This is because the only differ- 
ence with Vtrue is the gaussian error in Ar/ = .05 that affects 
v it! . The blue curves show this error is not a problem, be- 
cause the mode expansions are insensitive to gaussian noise 
in the 2500 galaxies, i.e. they are essentially perfect. This 
demonstrates that even though the TF noise is as large as 
for the actual data, the ability to find the correct flow, when 
characterized by only 20 numbers, is intact. 

Note also that the auto-covariance of (v tr ue — v ) (dot- 
dashed curves) is also is greatly reduced from that of i> itf . 
Recall that v & assumes linear theory estimated from the dis- 
tribution of 20000 galaxies. Occasionally the correlations 
are badly mistaken, when a large cluster (much larger than 
Virgo) is in the foreground and complicates the difference 
between physical and redshift space separations, but v is 
always an excellent approximation to the TF velocity. 
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Figure 9. Top: The velocity correlation of the real data and 15 
mock catalogs. The dashed red and solid red curves curve are, 
respectively, the correlations of v itl and v ia — v in the real data. 
The blue lines are each correlations of v itt — v for the mock 
catalogs. Bottom: Velocity correlations for 15 mock catalogs. The 
red curves are the velocity of Vi t f, the dot-dashed curves show 
the correlation of (vtrue — v g ), and the blue curves correspond 
to vtme — v it f). Both Vtrue and v g are first smoothed with the 
20 mode expansion before the autocovariance is computed. Note 
that the correlation of v itf — v is only slightly worse than the 
correlation of vtrue — v gs , showing that the velocity reconstruction 
dominates the errors. Note also that we are plotting the square 
root of the velocity correlation ^. 

6 THE CONSTRAINTS ON a, f3 

Equipped with the error covariance matrices, we proceed to 
minimize \ 2 m equation (|27p with respect to a and f). We 
shall present detailed results for a and /3 for fields expanded 
in 20 modes and 30 modes. 

The minimization is done by computing \ 2 on a grid of 
values in the plane a and j3. At the minimum point \ 2 = 
Xmin = 21.5 which is very reasonable given that the standard 
deviation from the expected value of 22 (20 mod e coefficients 
plus TF slope and zero point) is ~ \/H w 7 (|Press et al.l 
1 19921 ). 
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Figure 10. The difference A\ 2 = X 2 — Xmin versus P computed 
for three values of a, as indicated in the figure. Horizontal lines 
mark the 68% and 95.4 % CLs. 



Figure [10] shows the difference A\ 2 — \ 2 ~ Xmm as a 
function of f3 for three fixed values of a, as indicated in the 
figure. The horizontal lines indicate 68% and 95% confidence 
levels (CLs). Figure [12] is a contour plot of A\ 2 in the plane 
of a and j3. 

The higher frequency basis functions should probe 
smaller scales. Hence, if our assumption of linear bias which 
is independent of scale is valid then varying the number of 
modes, j m , should yield consistent constraints on f3. Figure 
llll shows the best fit j3 (thick red curve) and the correspond- 
ing 95% (2cr) CLs (thin blue solid lines) versus the number 
modes in the expansion, for a = 0. The circles show the 
"differential" best fit /3 obtained from a single mode as a 
function of the order of the mode. The 2a errorbars on this 
differential /3 are significantly enhanced beyond the second 
order mode. A few points lie at the ends of the errorbars 
corresponding to best fit /3 obtained at either 0.05 or 0.7 
which are the bounds of the range of /3 values used in the 
2MRS reconstruction. There is a hint that ft declines with 
increasing j m but this is completely dominated by the noise. 
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Figure 11. The thick solid red line is best fit /3 as a function of 
the number of modes in the velocity model, for a = 0. The thin 
blue lines mark the 95% CL on the best fit. The circles denote 
the "differential" best fit obtained with the single j th mode. 
After the first 4 modes, only 1 in 4 modes are represented. The 
errorbars attached to the circles are correspond to 95% CLs on 
the differential best fit. 
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Figure 12. Contour plot of A\ 2 in the plane of a and 0. The 
contours are 2.3, 6.17 and 9.2 corresponding to CLs of 68% 95.4% 
and 99%, respectively. 



This figure shows clearly that we get consistent constraints 
on f3 (the red curve) when varying the number of modes in 
the expansion. Further, it shows that most of the signal is 
contained the very few first modes. 



7 DISCUSSION 

The analysis reported demonstrates a good match between 
the 2MRS predicted and TF observed velocities. The analy- 
sis is unique in several respects. Firstly, it completely avoids 
dealing with covariance matrices of errors in the velocities 
estimated from the TF sample, Secondly, it uses elaborate 
mock galaxy catalogues to compute the error covariance in 
the predicted velocities from the redshift survey. Thirdly, 
the TF and predicted velocities are filtered in a very sim- 
ilar fashion, taking special care to minimize the effects of 
noise in the comparison. In this analysis faint galaxies with 
M > —20 are excised from the TF catalog since they system- 
atically show strong deviations from the linear TF relation. 

The comparison yields /3 — f/b = 0.33T0.04 (la error). 
The quoted error is not actually far from the limit of what 
the current data can constrain in the absence of any errors 
on the 2MRS predicted velocities 0. A moderate reduction 
of the errors by a factor of two requires a significant en- 
largement in the number of peculiar velocity measurements 
by a factor of 4, which could be done if dedicated time is 
available, but the TF samples already use the best local 
galaxies. Going to larger distance is not the answer, as the 
error of a peculiar velocity increases linearly with the dis- 
tance, and the 2MRS density field becomes very dilutely 
sampled. Two surveys, WALLABY, to be undertaken by 
the ASKAP telescope in Western Australia, and ALFALFA, 
an ongoing project at Arecibo, will hopefully produce good 
TF data for cz < 12000 km/s. Another strategy would in- 
volve peculiar velocities inferred from more precise distance 
indicators than the TF relation. SNe and surface bright- 
ness fluctuations techniques are likely candidates, but such 
measurements are available for a much smaller number of 
galaxies. Larger samples of local SNe are turning out to 
have increased errors (Ganeshalingam, private communica- 
tions), significan tly larger than previously measured (e.g. 
iRiess et al.1l 19971 ). 



2 In the absence of gravity errors and for close to the best fit 
value 0„ , the \ 2 function is approximated as 



X 



i 



B(0)al o f/al 



O and B = 

(/9//3 )(l + 1.5/3 )/(l + 1.5/3) approximates the dependence of 
v g on (in contrast to the dependence 0/0 o in reconstruction 
from galaxy distribution in real space). The lcr error in B is 

M s0/ ~ 0.06, where we took o- a = 0^/7 = 0.43 and 

X/j( a go) 2 = 52 as given from the solution with ft = 0.35. This 
error in B translates into an error of 0.03 in which is close to 
the error obtained with the full analysis. 
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The good match between the gravity and velocity 
fields implies that they probe the same underlying potential 
field within the framework of the gravitational instability 
paradigm for structure formation. The agreement is achieved 
assuming a linear biasing relation between mass and galaxies 
on large scales. Linear biasing is consist ent with theoretical 
predictions (e.g. iKauffmann et al.llT997l ) for the large scale 
clustering of galaxies. Further, no scale-dependent biasing 
seems to be required by the velocity comparison. However, 
the theo retically expected scale dependence of the bias f ac- 
tor (e.g. iPeacock fc Smithll2000l ; iDesiacques et al]|201Ch is 
well below the level which can be probed by the velocity 
comparison carried out here. There is also no clear indica- 
tion for a scale dependence bias fro m the observed ga laxy 
clustering on the relevant scales (e.g. IVerde et alJl2002h . 

DNW96, which compared the predicted velocities from 
the IRAS 1.2 Jy redshift surveys and the MARK III Tully- 
Fisher data, revealed systematic discrepancies that could not 
be attributed to errors in the data and the reconstruction 
methods. Inspection of the flow fields obtained in the current 
work with the those presented in DNW96 (see their figures 
9-13) clearly show that the problem lies in the MARK III 
data set. The velocity fields predicted from the IRAS and 
2MRS surveys have similar patterns which grossly deviate 
from MARK III but are in accordance with SFI++. The 
2MRS has all the attributes that one would want for estimat- 
ing the gravity field including a very weak bias. The survey 
was done by an instrument that was photometrically stable, 
which is important to avoid large scale drifts in the derived 
gravity field. The survey is far superior to the IRAS survey, 
the first full sky galaxy survey, which detected galaxies at 
60/im, a sign of star formation and not a good indicator of 
mass. However, that survey led to sensible results, and was 
not at fault for the disagreement 15 years ago. 

iRadburn-Smith et all (|2004h compared the predicted 
velocities from the IRAS PSCZ survey have also been com- 
pared with measured peculiar velocities of SNe. They found 
a best-fitting fi 1RAS = f/b IRAS = 0.55 ± 0.06. The lower 
value of (3 derived in our work could be due to a difference 
in the biasing factor between 2MRS and the IRAS galaxies, 
but we also emphasize that our estimation of the error in the 
predicted velocities sho uld be more reliable as it is based on 
realistic mock catalogs. iPike fc Hudson! l|2005h performed a 
comparison of the 2MRS predicted velocities with direct ve- 
locity measurements from three different samples, including 
836 SFI++ galaxies within cz = 5000 kms -1 . Their analysis 
yields f3 = 0.55 ± .05 for the comparison of gravity with the 
SFI++. They derive ft = 0.55 ± 0.05, inconsistent with our 
result at more than the 2.5a level. However, they did not 
calibrate their methods with advanced mock catalogues nor 
included the expected covariance of the predicted velocities. 

lLavaux et al.l (|2010h employed a sophisticated version 



of the nonlinear MAK reconstruction method l|Frisch et al.l 
l2002bh to compare the 2MRS predict ed velocities with the 
3K velocity catalog (|Tullv et alj|2008h of 1791 galaxies with 
redshifts < 3000 kms" 1 . They derive ft = 0.31 ±0.05, corre- 
sponding to f3 ~ 0.52. However, as they point out their error 
analysis is incomplete. Their method is promising as it takes 
intro account of nonlinear effects. Nevertheless, they do not 
account for the covariance of the errors in their smoothed 
observed velocities and predicted velocities. 

Both IPike fc Hud"sor] (|2005rj and lLavaux et ail (|2010h 
use iterative schemes based on lYahil et al.l 1 199l[ ) for de- 
riving the peculiar velocities from redshift surveys. These 
schemes rely on a relation between the peculiar velocity and 
density in real space. At any iteration, this relation is solved 
for new peculiar velocity given real space coordinates ob- 
tained from the observed redshifts by subtracting the old pe- 
culiar velocities derived in the previous iteration. We caution 
here that the these schemes are intrinsically biased: error in 
velocities used to estimate the distances will yield a biased 
density field in real space (see the [jX]for details). Hence the 
estimation of the velocity field is actually done from a biased 
distribution in real space. The bias produces an undesired 
smoothing of density field along the radial direction. The 
smoothing width (in kms -1 ) is equal to the rms random 
error in the velocities ~ 200 — 300 kms -1 feranchin i et alj 
l2002l;lNusser fc BranchinfeOOOl ') . Therefore, the bias is more 
pronounced in nonlinear methods which aim at probing 
small scales. 

Checks to find the best way of estimating the gravity 
field did not lead to improvements. Weighting the galaxy 
maps by the 2MRS luminosity led to a worse agreement, 
and recall that 2MRS is selected in K band, which is clos- 
est to a measure of the stellar mass. Giving the ellipti- 
cal galaxies double weigh t, as indicated by lensing analysis 
l|Mandelbaum et alj|200q ). did not improve the agreement. 
This in itself is not too surprising, because on large scale 
the 2MRS survey is dominated by spiral galaxies. It seems 
that the galaxies brighter than M* + 2 are each surrounded 
by a dark matter halo that has the same mass on average. 
There is no hint that the dark matter mass is larger if the 
luminosity is increased. 

Using our estimate for ft we can constrain the ampli- 
tude of mass fluctuations. As a measure of the amplitude 
we consider the rms of density fluctuations in spheres of 
8/i -1 Mpc in radius, denoted by ag and ag g for the mass 
and galaxy distribu tions, respectively. Adopting ft = 0.266 
llLarson et al.ll2010h gives /(ft, A = 1 - ft) = 0.483 ijLinderi 
l2005h Comparing this to our result = f/b = 0.33 ± 0.04 
(la error) we get a bias factor b — 1.46 ± .20 between 
the dark matter a nd the 2MRS ga l axy distribution. Taking 
a 8g = 0.97 ± 0.05 (|Westoved 120071 ; iReid et al.ll2010h . yields 
ag = agg/b = 0.65 ± 0.11 for the underlying mass density 
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field, marginally co nsistent with the latest WMAP results 
1 Larson et al.ll20ld ) of as = 0.8±0.03 (see also lJarosik et al.l 
( 20ld )). 



8 CONCLUSIONS 

We summarize the major conclusions of our work: 

• After a detailed examination of the 2MRS and SFI++ 
catalogs, we find the local gravity field to be a fine predictor 
of the local velocity field. Such a conclusion is a comfort 
for linear perturbation theory in an expanding universe and 
was certainly expected. It is interesting that the counts of 
galaxies give the best possible gravity field, reinforcing the 
old idea that the mass of the halo around a galaxy is not 
very well correlated with the luminosity of that galaxy. 

• We see no evidence that the dark matter does not follow 
the galaxy distribution, and is consistent with constant bias 
on large scales. There is no evidence for a non-linear bias in 
the local flows. A smooth component to the universe is not 
something testable with these methods. 

• Linear perturbation theory appears to be adequate for 
the large scales tested by our method. 

• The solution favors a = 0., no correlation between lu- 
minosity and mass, and /3 = . 33 ± .04, which is con sistent, 
but more than twice as tight aslErdogdu et al. | (|006al ). Using 
the derived Q from WMAP jjarosik et alj|20ld l. leads to an 
estimate erg = 0.65 ±0.10, deviant from WMAP's reports at 
the 1.5a level. 

• Our estimate of as gives the most precise value at z ~ 
and is useful for tests of the growth rate and Dark Energy. 

• The velocity-gravity comparison measures the acceler- 
ation on scales up to 30 — 50 Mpc. and since we derived a 
similar value of /3 as for clusters of galaxies, we conclude that 
dark matter appears to fully participate in the clustering on 
scales of a few Megaparsecs and larger. 

• We find no evidence for la rge-scale flow s such 
as reported by, for e xample jHudson et al.l Il999l : 
iHudson fc Ebelind Il997l ; iFeldman et alj l201Ch . Note 
that our analysis has not used the CMBR dipole, but we see 
a velocity field that is fully consi s tent w i th those previously 
repor ted (|Erdogdu et all 120061 . l006al ; lErdogdu fc Lahavl 
120091 ). which are consistent with the CMBR dipole radi- 
ation. We see no evidence that the dipole in the CMBR 
is produced by anything other than our motion in the 
universe. 
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APPENDIX A: A NOTE ON ITERATIVE 
SCHEME 

We caution of a possible systematic bias which may be 
important in iterative schemes for reconstructing velocities 
from the distribution of galaxies in redshift space. These 
schemes rely on the availability of a relation between the 
peculiar velocity and density in real space, e.g., the linear 
relation /divv = —5. At the end of any iteration interme- 
diate peculiar velocities are provided, which are used as to 
derive the distances from the redshifts in the next iteration. 
Given those distances, the adopted real space relation be- 
tween S and v is then solved to obtain a new guess for the 
velocities. The loop is continued until the change in the pe- 
culiar velocity between successive iterations becomes smaller 
than a certain threshold. Vanishing peculiar velocities could 
be taken as input for the first iteration. 

We demonstrate here that the real space distribution 
of galaxies as obtained from the output from any iteration 
scheme is biased. Hence the corresponding peculiar velocity 
is also biased. We will first show that a biased distribution 
in real space is obtained even if unbiased but noisy peculiar 
velocities are used to get the distances. We write the density 
of galaxies in real space at distance r in a given direction on 
the sky as 

n c (r)r 2 = (2na 2 )^ 1 ^ 2 J dss 2 n s e ^ ; (Al) 

where s = r+v is the radial redshift space coordinate, v(s) is 
the peculiar velocity of a galaxy present at s, and a is the rms 
of the error in the determination of the v(s). In the above we 
assume normal error distribution and a one-to-one mapping 
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between s and v, i.e. we neglect fingers-of-god effects and 
triple value zones. Working with the variable n = s — v(s) 
we get 

n c (r)r 2 = (27rcr 2 )~ 1/2 / dnr 2 n t (ri)e , (A2) 



where 
m(ri,) 



as 
dri 



(«(n)) , (A3) 



is the actual real space density at n. Therefore, errors in 
the peculiar velocities (even if unbiased relative to the true 
ones) cause a smearing of structure in the radial direction. 
This anisotropic smearing is important for scales ~ a (in 
kms -1 ). The bias is similar to the traditional inhomoge- 
neous Malmquist bias whi ch is usually encountere d in stud- 
ies of distance indicators (|Lvnden-Bell et al.lll98SI ). 

A self-consistent treatment of the bias should take into 
account the fact that v(s) used in lA2l and IA3I is the biased 
peculiar velocity obtained from n c . This could be done in the 
far observer limit, \v/r\ « 1, and for small perturbations 
where equation I A3I reduces to 

di> , . 

<5t = <5 S + — . (A4 
ds 



Substituting this into lA2l and Fourier transforming the result 
we get 

5 e (k) = [6,(k) + ik r v(k)] e~ k > 2/2 , (A5) 

where k r is the component of of k parallel to the line of sight 
and the tilde denotes quantities in k-space. Using this last 
equation in the linear 5 — v relation, fv = —i(k r /k 2 )5 e , we 
find 

* c (k) = d -^t (A6) 

V ' l + /(fc r /fc) 2 e- fe ^ 2 / 2 V ' 

instead of the usual unbiased expression obtained with a = 
0. 

A more complete analysis of the bias must incorporate 
the covariance of the errors in the derived peculiar velocity 
field. 



